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ABSTRACT 

The IceCube neutrino observatory pursues a follow-up program selecting interesting neutrino events in real¬ 
time and issuing alerts for electromagnetic follow-up observations. In March 2012, the most significant neu¬ 
trino alert during the first three years of operation was issued by IceCube. In the follow-up observations 
performed by the Palomar Transient Factory (PTF), a Type Iln supernova (SN) PTF12csy was found 0.2° away 
from the neutrino alert direction, with an error radius of 0.54°. It has a redshift of z = 0.0684, corresponding 
to a luminosity distance of about 300 Mpc and the Pan-STARRS 1 survey shows that its explosion time was 
at least 158 days (in host galaxy rest frame) before the neutrino alert, so that a causal connection is unlikely. 
The a posteriori significance of the chance detection of both the neutrinos and the SN at any epoch is 2.2 a 
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within IceCube’s 2011/12 data acquisition season. Also, a complementary neutrino analysis reveals no long¬ 
term signal over the course of one year. Therefore, we consider the SN detection coincidental and the neutrinos 
uncorrelated to the SN. However, the SN is unusual and interesting by itself: It is luminous and energetic, 
bearing strong resemblance to the SN Iln 2010jl, and shows signs of interaction of the SN ejecta with a dense 
circumstellar medium. High-energy neutrino emission is expected in models of diffusive shock acceleration, 
but at a low, non-detectable level for this specific SN. In this paper, we describe the SN PTF12csy and present 
both the neutrino and electromagnetic data, as well as their analysis. 

Keywords: circumstellar matter — galaxies: dwarf — neutrinos — shock waves — supernovae: individual 
(PTF12csy, SN 2010jl) 


1. INTRODUCTION 

IceCube is a cubic-kilometer-sized neutrino detector in¬ 
stalled in the ice at the geographic South Pole between depths 
of 1450 m and 2450 m (Achterberg et al. 2006). It consists 
of an array of 5160 photon sensors, called Digital Optical 
Modules (DOMs), attached to 86 cables, called strings. De¬ 
tector construction started in 2005 and finished in Decem¬ 
ber 2010. Neutrino observation relies on the optical detec¬ 
tion of Cherenkov radiation emitted by secondary particles 
produced in neutrino interactions in ice or bedrock near Ice- 
Cube. Due to the small neutrino interaction cross-section, 
the kilometer-scale detector has an effective area of only 
0.25 to 10m 2 for muon neutrinos of 1 to lOTeV energy 
(Aartsen et al. 2014b). As part of the Optical Follow-Up 
(OFU) program, the IceCube neutrino observatory records 
high-energy (~ lOOGeV... 1 PeV) neutrino events at a rate 
of about 3 mHz, about 250 per day. Those events are mostly 
(~ 90%) cosmic-ray induced neutrinos from the atmosphere, 
referred to as atmospheric neutrinos, with about a 10% con¬ 
tamination of cosmic-ray induced atmospheric muons. 

Routine neutrino analyses in IceCube, which are referred 
to as offline analyses, are performed after a certain amount 
of data, e.g. one or several years, has been collected. They 
benefit from events being put through computationally ex¬ 
pensive reconstructions, as well as information on detector 
performance that become available only days or weeks after 
data acquisition. In contrast, neutrino analyses running on¬ 
line will not have access to such information, but have the 
advantage of being near real-time—results are available with 
a latency of ~ 3 min. With such a short latency neutrino analy¬ 
sis, multi-wavelength follow-up observations can be triggered 
by neutrino events. These follow-up data have the potential 
to reveal the electromagnetic counterpart of a transient neu¬ 
trino source, which might otherwise be missed and thus be 
unavailable for further observations. In addition, the coinci¬ 
dent detection of neutrino and electromagnetic emission can 
be statistically more significant and provide more informa¬ 
tion about the physics of the source than the neutrino detec¬ 
tion alone. Another advantage of an online analysis is the 
prompt availability of the reconstructed neutrino dataset and 
thus the possibility of fast response analyses. Thus, IceCube’s 
online neutrino analysis efforts have also enabled fast y-ray 
burst (GRB) searches like the one following GRB 130427A, 
published in a GCN Circular (Blaufuss 2013). 

The online search for short transient neutrino sources (or¬ 
der of < 100 s) is mostly motivated by models of neutrinos 
from long duration GRBs (Waxman & Bahcall 1997; Murase 
2008; Murase et al. 2006; Becker et al. 2006) and from choked 
jet supernovae (SNe) (Razzaque et al. 2004; Ando & Bea- 
com 2005). The two source classes are related: Both are 

^ Corresponding author (voge@physik.uni-bonn.de). 


thought to host a jet, which is highly relativistic in case of 
long GRBs, but only mildly relativistic in case of choked jet 
SNe. Long GRB progenitors are conceived to be Wolf-Rayet 
stars that have lost their outer hydrogen and helium envelope 
(Meszaros 2006), while choked jet SNe can still have those 
outer layers (Ando & Beacom 2005) which are important for 
efficient high-energy (HE) neutrino production. The choked 
jet is more baryon-rich and has a much lower Lorentz factor 
r«3 than the GRB jet with T > 100. It cannot penetrate the 
stellar envelope and remains optically thick, making it invis¬ 
ible in y-rays. The neutrinos produced at TeV energies can 
escape nevertheless and may trigger the discovery of the SN 
in other channels. In particular, for a bright nearby SN, a neu¬ 
trino detection would enable the acquisition of optical data 
during the rise of the light curve, strengthening the time cor¬ 
relation between the neutrino burst and the optical supernova 
via the improved estimate of the explosion time (Cowen et al. 
2010). Mildly relativistic jets may occur in a much larger 
fraction of core-collapse SNe than highly relativistic jets, i.e. 
GRBs (Razzaque et al. 2004; Ando & Beacom 2005). 

Detections of neutrinos from GRBs and choked jet SNe 
would be a remarkable discovery and could provide impor¬ 
tant insight into the supernova-GRB connection and the un¬ 
derlying jet physics (Ando & Beacom 2005). Both sources 
are expected to emit a short, about 10 s long burst of neutrinos 
(Ando & Beacom 2005) either 10 to 100 s before or at the time 
of the y-ray burst (if detectable) (Meszaros 2006), setting the 
natural time scale of the neutrino search. After recording the 
neutrino burst, follow-up observations can be used to iden¬ 
tify the counterpart of the transient neutrino source. A y-ray 
burst can be identified either via the prompt y-ray emission 
lasting up to about 150s (Baret et al. 2011) or via the opti¬ 
cal and X-ray afterglow lasting up to several hours (Gehrels 
et al. 2004). The latter involves instruments of limited field of 
view and thus requires telescope slew, but has the advantage 
of much better angular resolution well below an arc minute 
(Gehrels et al. 2004). A fast response within minutes to hours 
is required for a GRB afterglow follow-up. A choked jet SN 
is found by detecting a shock breakout or a SN light curve in 
the follow-up images, slowly rising and then declining within 
weeks after the neutrino burst. Following this scientific moti¬ 
vation, an online neutrino analysis, targeted at SN and GRB 
afterglow detection, was installed at IceCube in 2008 (see 
Section 2). 

In addition to the transient neutrino emission within 100s, 
as discussed above, SNe can be promising sources of high- 
energy neutrinos over longer time scales. In this paper, the 
class of Type Iln SNe (Schlegel 1990; Filippenko 1997) is 
explored further. These are core-collapse SNe embedded in a 
dense circumstellar medium (CSM) that was ejected in a pre¬ 
explosion phase. Following the explosion, the SN ejecta plow 
through the dense CSM and collisionless shocks can form and 
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accelerate particles, which may create high-energy neutrinos. 
This is comparable to a SN remnant, but on a much shorter 
time scale of 1 to lOmonths (Murase et al. 2011; Katz et al. 
2012 ). 

SNe Iln (“n” for narrow) are spectrally characterized by the 
presence of strong emission lines, most notably Ha, that have 
a narrow component, together with blue continuum emission 
(Schlegel 1990; Filippenko 1997). The narrow component is 
interpreted to originate from surrounding Hll regions, and 
the generally slow spectral evolution is due to the presence 
of a high-density CSM (Schlegel 1990). Since interaction 
of the SN ejecta with the dense CSM can lead to the con¬ 
version of a large fraction of the ejecta’s kinetic energy to 
radiation, SNe Iln are on average more luminous than other 
SNe II (Richardson et al. 2002). They generally fade quite 
slowly and some belong to the most luminous SNe (Filip¬ 
penko 1997). There is much diversity within the subclass of 
SNe Iln (Filippenko 1997; Richardson et al. 2002), both spec¬ 
troscopically and photometrically, which can be explained by 
a diversity of progenitor stars and mass loss histories prior to 
explosion (Moriya & Tominaga 2012). Recently, there have 
been observations of eruptions prior to SN Iln explosions as¬ 
sociated with mass loss which explain the existence of the 
dense CSM shells (e.g. Ofek et al. 2014a). 

In this paper, we report the discovery of a Type Iln SN 
in the optical observations triggered by an IceCube neutrino 
alert from March 2012, and analyze the available neutrino 
and electromagnetic observations. The SN is already at a late 
stage at the time of the neutrino detection, which means that 
the neutrino-SN connection is presumably coincidental. The 
paper is organized as follows: Section 2 introduces the op¬ 
tical follow-up system of IceCube. Section 3 gives details 
about the neutrino alert triggering the follow-up observations 
that led to the SN discovery. Section 4 reports limits from a 
complementary offline neutrino search and X-ray limits. The 
UV and optical data that were obtained are discussed in depth 
in Section 5. We finally summarize the results and give a con¬ 
clusion in Section 6. 

2. THE OPTICAL AND X-RAY FOLLOW-UP SYSTEM 

In late 2008, an online neutrino event selection was set up at 
IceCube, looking for muons produced by charged current in¬ 
teractions of muon neutrinos in or near the IceCube detector. 
The analysis is running in real-time within the limited com¬ 
puting resources at the South Pole, capable of reconstructing 
and filtering the neutrinos and sending alerts to follow-up in¬ 
struments with a latency of only a few minutes (Aartsen et al. 
2013; Abbasi et al. 2012; Kowalski & Mohr 2007). 

The optical (OFU) and X-ray (XFU) real-time follow-up 
programs currently encompass three follow-up instruments: 
the Robotic Optical Transient Search Experiment (ROTSE) 
(Akerlof et al. 2003), the Palomar Transient Factory (PTF) 
(Law et al. 2009; Rau et al. 2009) and the Swift satellite 
(Gehrels et al. 2004). These triggered observations were 
supplemented with a retrospective search through the Pan- 
STARRS1 3 71 survey data (Kaiser 2004; Magnier et al. 2013), 
which is discussed further in Section 5. In addition, there 
is also a real-time y-ray follow-up program (GFU) targeting 
slower transients (time scale of weeks), e.g. flaring Active 
Galactic Nuclei, that is sending alerts to the y-ray telescopes 
MAGIC and VERITAS (Aartsen et al. 2013). 

The background of cosmic-ray induced muons from the at¬ 
mosphere above the detector amounts to ~ 10 6 muon events 
per neutrino event. In a first step, it is reduced by limiting 


the sample to the northern hemisphere, using the Earth as a 
muon shield and selecting only muon tracks that are recon¬ 
structed as up-going in the detector. Afterwards, cuts on qual¬ 
ity parameters similar to those in Aartsen et al. (2014b) are 
applied to reject mis-reconstructed muon events: After an ini¬ 
tial removal of likely noise signals, a chain of reconstructions 
is performed, which utilizes the spatial and temporal distribu¬ 
tion of recorded photons on the DOMs. Starting with a simple 
linear track algorithm, more advanced reconstructions are em¬ 
ployed, seeded with the respective preceding reconstruction. 
The advanced reconstructions maximize a likelihood that ac¬ 
counts for the optical properties of the ice for photon propaga¬ 
tion (Ahrens et al. 2004). The final reconstruction in this chain 
is used to select events that are up-going in the detector. At 
this point, because of the vast amount of atmospheric muons, 
the data are still dominated by (down-going) muons that are 
mis-reconstructed as up-going. To reject this remaining back¬ 
ground, high quality events are selected, where the selection 
parameters are derived from the value of the maximized track 
likelihood and from the number and geometry of recorded sig¬ 
nals with a detection time compatible with unscattered photon 
propagation. Alternatively, events with a large number of to¬ 
tal recorded photon signals are selected. The resulting analy¬ 
sis sample has an event rate of about 3 mHz, of which ~ 90% 
are atmospheric muon neutrinos that have passed through the 
Earth and ~ 10% are contamination of cosmic-ray induced 
muons from the atmosphere. 

In the analysis sample, the median angular resolution of 
the neutrino direction is about 1° for a multi-TeV muon neu¬ 
trino charged current interaction event, and 0.6° or less for 
lOOTeV and higher energies. The angular resolution of the 
sample is estimated using Monte Carlo (MC) simulation. Ad¬ 
ditionally, an estimator of the directional uncertainty is com¬ 
puted for each event, which is based on the shape of the 
reconstruction likelihood close to the found maximum (Ne- 
unhoffer 2006). It is calibrated such that its median matches 
the MC derived angular resolution. We define the bulk of a 
sample as events within the central 90% of the energy dis¬ 
tribution. The bulk of the main background contained in the 
analysis sample, atmospheric neutrinos, has energies between 
160GeV and 7TeV. In contrast, the bulk of signal events 
from an unbroken E~ 2 power-law spectrum would have en¬ 
ergies 1.2TeV < E v < 1.2PeV. Signal neutrinos from GRBs 
are expected to follow a spectrum similar to E ~ 2 , which has 
cut-off energies between ~ 1 PeV to ~ lEeV (Murase & Na- 
gataki 2006). Choked jet SNe are predicted to have lower cut¬ 
off energies around 20TeV (Ando & Beacom 2005) so that 
the overlap with atmospheric neutrinos is presumably much 
larger, while Iln neutrinos may have higher cut-offs of 70 to 
200TeV (see Section 4.1). 

In order to suppress background from atmospheric neutri¬ 
nos, a multiplet of at least two neutrinos within 100 seconds 
and angular separation of 3.5° or less is required to trigger an 
alert. In addition, since 2011 mid-September, a test statistic 
is used, providing a single parameter for selection of the most 
significant alerts. It was derived as the analytic maximiza¬ 
tion of a likelihood ratio following Braun et al. (2010), for the 
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special case of a neutrino doublet with rich signal content: 


Avp 2 

A = —— +2]n(2n(jq) 

°q 


— 21n 


^ 1 — exp 



+ 2 In 



( 1 ) 


where the time between the neutrinos in the doublet is de¬ 
noted as AT, and their angular separation as A 1 ?. The quanti¬ 
ties cr 2 = of + of and of = (l/of + 1/ of) 1 depend on the 
event-by-event directional uncertainties Oi and Oh of the two 
neutrino events, typically ~ 1°. The angle Q,\ corresponds to 
the circularized angular radius of the field of view (FoV) of 
the follow-up telescope. It is set to 0.5° for Swift and 0.9° for 
ROTSE and PTF. 

The test statistic A is smaller for more signal-like alerts, 
which have small separation ATk small time difference AT 
and a high chance to lie in the FoV of the telescope. Thus, A 
is a useful parameter to separate signal and background alerts. 
For each follow-up program, a specific cut on A is applied 
in order to send the most significant alerts to the follow-up 
instruments. In the 2011/12 data acquisition (DAQ) season, 
which is discussed here, a cut of A < —7.4 was used for the 
ROTSE follow-up, while cut values of —10.3 and —8.8 were 
used for PTF and Swift. Multiplets of multiplicity higher than 
two are passed directly to all follow-up instruments. Since 
the expected background rate is low (~ 0.03 per year), each 
observation of a triplet or higher order multiplet is significant 
by itself. 

ROTSE (Akerlof et al. 2003) is a network of four optical 
telescopes with 0.45 m aperture and 1.85° x 1.85° FoV, lo¬ 
cated in Australia, Texas, Namibia and Turkey. Since late 
2012, only the two northern hemisphere telescopes continue 
operation. ROTSE is a completely automatic and autonomous 
system that can receive alerts, perform observations and send 
resulting data without requiring human interaction. The lim¬ 
iting magnitude of ~ 16 to 17 mag is however insufficient 
to discover faint or far SNe. For instance, for a very bright 
SN with —20 mag absolute magnitude, the detection radius is 
about 160 to 250Mpc, while a faint SN with — 17mag is only 
visible within a radius of 40 to 65Mpc. IceCube has been 
sending ~ 25 alerts per year to ROTSE, since 2008 Decem¬ 
ber. The first 116 alerts, with a background expectation of 
104.7 ± 10.2 alerts, were followed up with a median latency 
of 27.2 hours between the neutrino alert and start of the first 
follow-up observation. 

PTF (Law et al. 2009; Rau et al. 2009) is a survey based 
at the Palomar Observatory in California, USA. It utilizes the 
1.2 m Oschin Schmidt telescope on Mount Palomar. The fo¬ 
cal plane is equipped with a mosaic of 11 CCDs with field of 
view of 7.26deg 2 . The typical //-band limiting magnitude of 
PTF during dark time is about 21 mag. All the PTF data are re¬ 
duced using the LBNL real-time pipeline responsible for tran¬ 
sient identification and the IPAC pipeline described in Laher 
et al. (2014). The image photometric calibration is described 
in Ofek et al. (2012). PTF pursues a number of science goals, 
most notably the discovery and observation of SNe. Several 
other telescopes in Palomar and at other locations can be used 
for photometric and spectroscopic follow-up observation. Ice- 
Cube has been sending ~ 7 alerts per year to PTF, since Au¬ 
gust 2010. The first 23 alerts, with a background expectation 
of 19.2 ±4.4 alerts, were followed up with a median latency 
of 34.9 hours between the neutrino alert and start of the first 


follow-up observation. 

Swift (Gehrels et al. 2004) is a satellite operated by NASA 
and boards various instruments: a 170 to 600nm ultravio¬ 
let/optical telescope (UVOT), a 0.3 to lOkeV X-ray telescope 
(XRT) and a 15 to 150keV hard X-ray Burst Alert Telescope 
(BAT). Swift’s main goal is the discovery and study of GRBs, 
of which it detects about 100 per year (Lien et al. 2014) (one 
third of all GRBs) 1 . IceCube’s X-ray follow-up program trig¬ 
gers Swift’s XRT, which can provide valuable information by 
observing a GRB afterglow in X-rays. The XRT has a FoV of 
only 0.4° in diameter, hence Swift performs seven pointings 
for each IceCube follow-up, resulting in an effective FoV of 
about 1° in diameter. IceCube has been sending ~ 6 alerts per 
year to Swift, since February 2011. The first 18 alerts, with a 
background expectation of 18.0±4.2 alerts, were followed up 
with a median latency of 1.9 hours between the neutrino alert 
and start of the first follow-up observation (see Evans et al. 
2015). 

3. NEUTRINO ALERT AND DISCOVERY OF PTF12CSY 

On 2012 March 30 (MJD 56016), the most significant 
alert since initiation of the follow-up program (significance 
of ~ 2.7 <7, converting the A cumulative distribution function 
value to single-sided Gaussian std. deviations) was recorded 
and sent to ROTSE and PTF simultaneously. The significance 
was also above the threshold for Swift (~ 1cr), however it was 
within Swift’s moon proximity constraint 2 , which delayed the 
observations by three weeks. The two neutrino events causing 
the alert happened on 2012 March 30 at 01:06:58 UT (MJD 
56016.046505) and 1.79 seconds later, with an angular sep¬ 
aration of 1.32°. The combined average neutrino direction 
is at right ascension 6 h 57 ra 45 s and declination 17°11’24” in 
J2000 with an error radius of a w = 0.54°. This average is a 
weighted arithmetic mean, weighting the individual directions 
with their inverse squared error, given by the event-by-event 
directional uncertainty (s.a.). The error a w is defined after 
Eq. 1. This assumes that the two neutrino events were emitted 
by a point source at a single fixed position. A variance of the 
individual true neutrino directions does not need to be taken 
into account, since the assumed intrinsic variance is zero in 
case of a point source. This leads to a relatively small error 
on the average direction. 

The main event properties are summarized in Table 1 : the 
occurrence time on 2012 March 30, the reconstructed muon 
energy proxy E^ (see Aartsen et al. 2014), and the estimated 
directional error a>p. The quantity is a fit parameter and 
serves as a proxy for the muon energy, however it is not an 
estimator of the true muon energy. The energy E v of the neu¬ 
trino that produced the muon is not directly observable, since 
only the muon crossing the detector is accessible. However, 
using Monte Carlo (MC) simulated neutrino events, one can 
use the muon energy proxy E f j to compare with MC events 
having a similar E^ value. From those MC events, a distri¬ 
bution of the true muon energy E v can be derived, which de¬ 
pends on the assumed underlying neutrino energy spectrum. 
E v probability density functions are filled with true energy 
values of MC events that have the reconstructed muon energy 
proxy E u not more than 10% and the reconstructed zenith an¬ 
gle cos(O) not more than 0.1 away from the observed values 

1 http://gcn.gsfc.nasa.gov/ and http://grbweb.icecube. 
wise.edu 

2 Swift is unable to observe sources closer than 15° to the moon. 
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Table 1 

Properties of the neutrino alert events 


Time 

(UT) 

Off) 

V 

(GeV) 

E v b (Atm.) 
(TeV) 

E v b ( E ~ 3 ) 
(TeV) 

E v b (E- 2 ) 
(TeV) 

01:06:58 

01:07:00 

0.96 

0.66 

1155 

3345 

U3 —0.4 

0 9+ 6 - 7 
uy -0.7 

0 7+ 5 - 6 
—0.5 

1 5+ 14 - 8 

A ’ —1.3 

c 4+292.0 

ic 7 +611.5 

1 ' > ' 7 —14.5 


a En is only a proxy correlated with muon energy, but not an estimator of the 
true muon energy. 

b E v is median neutrino energy with 90% C.L. error interval. 



Right Ascension [°] 


Figure 1. Map of the sky with the two neutrino event directions, the average 
neutrino direction, and the location of SN PTF12csy. Estimated reconstruc¬ 
tion errors are indicated with circles, the PTF FoV is shown as dashed box. 
The positions of the PTF survey camera CCD chips are plotted with dot¬ 
ted lines and the chip number is printed on each chip’s field (cf. Law et al. 
(2009)). Note that chip 03 is not operational and thus hatched in the plot. 

of the two alert events. The median and the central 90% C.L. 
interval are calculated from the E v probability density func¬ 
tions. The results are listed in Table 1, where the assumed 
neutrino spectrum is given in parentheses. 

Follow-up observations at the direction of the neutrino 
alert were performed with multiple instruments (see Sec¬ 
tion 5.1). In the PTF images, a core-collapse super¬ 
nova, named PTFUcsy, was discovered at right ascension 
6 h 58 m 32 s .744 and declination 17°15’44.37” (J2000), only 
0.2° away from the average neutrino direction, see Figs. 1 
and 2. This was a promising candidate for the source of the 
neutrinos, but a search of the Pan-STARRSl archive (see Sec¬ 
tion 5.1) revealed that it was already at least 169 observer 
frame days old, i.e. 158 days in host galaxy rest frame, at 
the time of the neutrino alert. Therefore, it is highly unlikely 
that the neutrinos were produced by a jet at the SN site, as this 
is expected to happen immediately after core collapse in the 
choked jet scenario (Ando & Beacom 2005). 

However, steady neutrino emission on a time scale of sev¬ 
eral months is a possibility and explored in Section 4.1 . 

3.1. Significance of Alert and SN Detection 

The value of the test statistic A for the neutrino doublet 
amounts to —18.1. The background distribution of A is 
constructed from experimental data, containing mostly at¬ 
mospheric neutrinos, by randomly permuting (shuffling) the 
event times and calculating equatorial coordinates, i.e. right 
ascension and declination, from local coordinates, i.e. zenith 
and azimuth angle, using the new times. That way, all detec- 


NEW REF 



Figure 2. New image, reference image and post-subtraction image of the 
PTF discovery of PTF12csy from 2012 April 09, with the location of 
PTF12csy in the center. This image shows only a small fraction of the PTF 
FoV. The image from the Sloan Digital Sky Survey (SDSS-III) DR 12 (Eisen- 
stein et al. 2011; Gunn et al. 2006; Ahn et al. 2014; Alam et al. 2015) is shown 
for reference, showing a faint host galaxy to the south of the SN. 


tor effects are entirely preserved, e.g. the distribution of the 
azimuth angle, which has more events at angles where de¬ 
tector strings are aligned, and the time distribution, which is 
affected by seasonal variations. At the same time, all potential 
correlations between the events in time and space, and thus a 
potential signal, are destroyed. 

The false alarm rate (FAR) for an alert with A < —18.1 
is 0.226 yr , calculated via integration of the A distribu¬ 
tion below —18.1. Considering the OFU live time of 220.1 
days in the data acquisition season of the alert, September 
2011 to May 2012, yields N (A < —18.1) = 0.136 false alerts. 
Hence, the probability, or p-value, for one or more alerts 
at least as signal-like to happen by chance in this period is 
1 — Pp o isson(0;(V(A < —18.1)) ss 12.7%. The OFU system had 
already been sending alerts to PTF for ~ 460 days at the time 
of the alert. Scaling up the number of expected alerts with 
A < —18.1, one derives a probability of ~ 24% during 460 
days. 

The estimated explosion time of SN PTF12csy does not fall 
within the a priori defined time window for a neutrino-SN 
coincidence of 0(1 day). It is thus not considered an a priori 
detection of the follow-up program. Despite this fact, for il¬ 
lustrative purposes, we calculate the a posteriori probability 
that a random core-collapse SN (CCSN) of any type, at any 
stage after explosion, is found coincidentally within the error 
radius of this neutrino doublet and within the luminosity dis¬ 
tance of PTF12csy, i.e. 300Mpc. The number of such random 
SN detections is 


-V det 


/-aooMpc dN SN 

An Jo dtdV 


T (ffijim ,M,r) Anr 2 dr 


( 2 ) 


where is the solid angle of the doublet error circle (blue 
circle in Fig. 1), which is ~ 0.93 (°) 2 . For the volumetric 
CCSN rate dN S x/(dtdV), a value of 0.78 x lO^Mpc -3 ^- 1 
is used, Horiuchi et al. (see 2013, sec. 4.1). The control time 
T (,wiimjM, r) is the average time window in which a SN is de- 
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tectable, i.e. brighter than the limiting magnitude. It depends 
on the distance to the source r, the peak absolute magnitude 
M of the SN, the limiting magnitude ni\\ m of the telescope, and 
the shape of the light curve which is adopted from a SN tem¬ 
plate web page by P. E. Nugent 3 . It is assumed that M follows 
a Gaussian distribution with mean —17.5 mag and standard 
deviation cr = 1 mag, based on Richardson et al. (2002). For 
PTF, mii m = (19.5 ± 1.0) mag is assumed. 

The resulting expectation value for coincidental SN detec¬ 
tions is Ndet ~ 0.016, which results in a Poisson probability 
of ~ 1.6% to detect any CCSN within the neutrino alert’s er¬ 
ror radius. Combining this probability with the probability of 
12.7% for the neutrino alert, Fisher’s method (Brown 1975; 
Littell & Folks 1971; Fisher 1950, 1990) delivers a combined 
p-value of 1.4%, corresponding to a significance of 2.2 (7. For 
the total live time of 460 days, the combined p-value is 2.5% 
which corresponds to a 2.0(7 significance. This means, even 
ignoring the a posteriori nature of the p-value, a chance coin¬ 
cidence of the neutrino doublet and the SN detection cannot 
be ruled out and thus we consider the SN detection to be co¬ 
incidental. 

The following section reports about the available high- 
energy follow-up data. Limits on a possible long-term neu¬ 
trino emission from PTF12csy are set using one year of Ice- 
Cube data. Limits on the X-ray flux were obtained using the 
Swift satellite. Section 5 deals with the analysis of the low- 
energy optical and UV data as the SN detection is significant 
and interesting by itself. 

4. HIGH-ENERGY FOLLOW-UP DATA 
4.1. Offline Analysis of Neutrino Data 

Type Iln supernovae, such as PTF12csy, are a promising 
class of high-energy transients (see Murase et al. (2011)). 
The expected duration of neutrino emission from SNe Iln is 
1 to 10 months, hence it is extremely unlikely that two neu¬ 
trinos arrive within less than 2 s, so late after the SN explo¬ 
sion. However, to test the possibility of a long-term emission, 
a search for neutrinos from PTF12csy within a search window 
of roughly one year is conducted. 

After the core-collapse of a SN Iln, the supernova ejecta 
are crashing into massive circumstellar medium (CSM) shells, 
producing a pair of shocks: a forward and a reverse shock. 
Cosmic rays (CRs) may be accelerated and multi-TeV neutri¬ 
nos produced, potentially detectable with IceCube. The colli¬ 
sionless shocks generating the neutrinos are expected to gen¬ 
erate X-rays as well at late times (see e.g. Katz et al. (2012); 
Svirski et al. (2012); Ofek et al. (2013)), but no X-rays were 
detected for PTF12csy, likely because of the large distance to 
the SN. 

Following Murase et al. (2011) and Murase et al. (2014) 
(see also Katz et al. 2012), we model high-energy neutrino 
emission from PTF12csy. As a simplified approach, in order 
to get an order-of-magnitude estimate of the expected event 
rate, we perform the following calculation: The CSM density 
profile is calculated using Murase et al. (2014, eq. 4, eq. A4). 
The proton spectrum is modeled with a power law index —2, 
as in Murase et al. (2014, eq. A7), with a cut-off energy given 
by the maximum energy of accelerated protons. The latter is 
determined by comparing the proton acceleration time scale 
either with the dynamical time scale (Murase et al. 2014, 
eq. 28), or, if pp energy losses are relevant, with the cooling 

3 https ://c3 .lbl.gov/nugent/nugent_templates.html 


time scale (Murase et al. 2014, eq. 30). The lower of the two 
maximum energies is the one that needs to be considered. The 
proton spectrum is normalized to the total CR energy £cr by 
assuming that a fraction of the kinetic energy of the ejecta E e j 
is converted into CRs, that is £cr = CcR^ej (compare Murase 
et al. 2014, eq. 3, eq. 25, using CSM shell mass Mqsm 
SN ejecta mass M e j). Other model parameters are the break¬ 
out radius A),,, and shock velocity v s h oc k- The expected neu¬ 
trino spectrum from pp interaction is derived from the semi- 
analytical description in Kelner et al. (2006), taking into ac¬ 
count the meson production efficiency (Murase et al. 2014, 
eq. 35, eq. 36). It is distance-scaled and folded with IceCube’s 
effective area from Aartsen et al. (2014b) to obtain the ex¬ 
pected number of events. It is found that inserting commonly 
assumed values (Murase et al. 2014; Margutti et al. 2014; 
Murase et al. 2011) of E e j = lOFboi = 2.1 x 10 51 erg (with 
the bolometric energy found in Section 5.2.5), £cr = 0.1, 
Mgj = 10M e , Rbo ~ 6 x 10 15 cm, and v s h 0 ck ~ 5000kms" 1 , on 
average only 0.07 IceCube neutrino detections are expected. 

Despite the low expectation for the neutrino fluence, we 
search for a long-term neutrino signal from PTF12csy in the 
IceCube data. As a more elaborate approach compared to the 
simplified approximation described above, we test the neu¬ 
trino emission models A and B given in Fig. 1 of Murase et al. 
(201 1), which are two representative cases of CR accelerating 
scenarios: Model A corresponds to a CSM shell with a high 
density of 1 x 10 n cm~ 3 at a small radius of 1 x 10 15 5 cm, 
while model B is the opposite with a density of 1 x 10 7 5 cm~ 3 
at radius 1 x 10 16 5 cm. For the ejecta, a kinetic energy of 
1 x 10 51 erg, a velocity of 1 x 10 4 kms _1 , and a mass of sev¬ 
eral solar masses, lower than the CSM mass, are assumed. 
Model A is close to a scenario explaining superluminous SNe 
Iln such as SN 2006gy, while model B is a good description 
for dimmer, but longer lasting SNe like SN 2008iy. Both mod¬ 
els have a neutrino energy spectrum close to E~ 2 , with a cut¬ 
off energy around 70TeV for model A, around 84TeV for the 
forward shock (FS) in model B, and around 275 TeV for the 
reverse shock (RS) in model B. In model A, only the reverse 
shock is of importance for CR acceleration. The suggested 
emission time scales are 1 x 10 7 s = 115 d for model A and 
1 x 10 7 - 8 s = 366d for model B (Murase et al. 201 1). 

Figure 3 shows the model fluences scaled down to a lumi¬ 
nosity distance of 308 Mpc. We analyze about one year of 
IceCube data: the entire IceCube 86 strings data acquisition 
season 2011/12, from 2011 May 13 to 2012 May 15. The 
long search window is motivated by the large uncertainty on 
the explosion date (between 2011 March 21 and 2011 Octo¬ 
ber 13) as well as the long duration of neutrino emission for 
some scenarios like ~ 700 days for model B in Murase et al. 
(2011). For simplicity, we assume that the entire fluence was 
emitted during the 1 yr search window. We use the neutrino 
sample of the IceCube optical follow-up system and perform 
a statistical point source analysis of neutrino events close to 
the position of the supernova, based on Braun et al. (2008). 

Each neutrino candidate event i is given both a signal and 
background probability .S', and If which are combined in the 
likelihood function 



The variable n s is the number of signal events contained in 
the sample, which is fitted to maximize the likelihood, and 
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Figure 3. Neutrino fluence at Earth from PTF12csy (solid lines) and derived 
upper limits set by IceCube (dashed lines, corresponding gray scales) as func¬ 
tion of energy for the tested models A, B reverse shock (RS), and B forward 
shock (FS) from Murase et al. (2011). 


N is the total number of selected neutrino candidate events. 
The signal probability Sj is the product of the spatial Gaus¬ 
sian probability density function (PDF) and the energy PDF 
P(Ej\<p), i.e. the probability of a signal event having recon¬ 
structed energy E, given the neutrino spectrum <j> of the source 
(derived from Monte Carlo simulation): 

Si(ai,Xi,x s ,Ei) = exp ( ^2 a 2^ ) p i E i\<l>)- ( 4 ) 


Here, cr, is the event’s angular error estimate, x,- is the event’s 
reconstructed direction, and x s the SN or neutrino source posi¬ 
tion. The background probability B, contains the energy PDF 
and a normalization constant for the background from atmo¬ 
spheric neutrinos. 

We define the test statistic as likelihood ratio (Braun et al. 
2008) 


A 


-2 log 


-g(O) ' 

■&{n s ) \ 1 


(5) 


serving as a powerful test for separating the null hypothesis 
from the hypothesis of signal event contribution. Here, h s is 
the number of signal events that maximizes the likelihood and 
corresponds to the most likely description of the data. 

The result of the maximum likelihood fit is n s = 0 both for 
model A and B, i.e. we see no sign of signal contribution in 
our neutrino event sample. We set 90% C.L. Neyman upper 
limits (see Neyman (1937), reprinted in Neyman (1967)) on 
the tested neutrino fluence models, which amount to ~ 1500 
and ~ 1300 times the fluences given for model A and B above, 
respectively. The limits are much higher than the fluence pre¬ 
diction because of IceCube being insensitive to SNe Iln at 
such large distances. Figure 3 shows a plot of the tested neu¬ 
trino fluence and the limits set using 1 yr of IceCube data. 

This null result and the large distance to the SN further sup¬ 
port our conclusion that the SN detection was coincidental. 

Nevertheless, it is interesting to roughly estimate the hypo¬ 
thetical emitted neutrino fluence: We take the median neu¬ 
trino energies of the two alert neutrinos from Table 1 and 
look up the effective areas for the OFU neutrino sample at 
the respective energies. From this, we derive a hypothetical 
neutrino fluence of 3.2 x 10 4 erg cm 2 for a source spectrum 
°c E - 2 (E v = 5.4TeV, 15.7TeV), or 10.8 x 10~ 4 ergcm- 2 for 
a source spectrum ( E v = 0.7 TeV, 1.5 TeV). Assum¬ 


Table 2 

Swift XRT observations of PTF12csy 


Time (MJD) 

Exposure (ks) 

c 

Cbg 

CUL 

‘I’UL 

56037.15 

4.9 

1 

1.47 

1.3 

4.6 

56245.29 

2.0 

1 

0.62 

2.1 

7.4 

56246.04 

1.2 

2 

0.44 

9.8 

30.0 

56247.62 

5.0 

2 

1.35 

2.1 

7.4 

Sum 

13.0 

6 

3.71 

1.3 

4.6 


Note. — Energy range: 0.2 to lOkeV. c: measured 
counts within a 28" aperture. r| n ,: expected background 
counts within the same aperture, cul: 3(7 upper limit on 
the X-ray count rate in 10 A 1 . upper limit on the 

unabsorbed source flux in 10 l4 crgcm " s 1 . 

ing that this neutrino fluence was emitted by the SN, this 
would imply a radiated neutrino energy of ~ 3.4 x 10 51 erg or 
~ 1.2 x 10 52 erg using the luminosity distance of ~ 300Mpc, 
corresponding to about 15 or about 50 times the radiated 
electromagnetic energy of £boi = 2.1 x 10 5 °erg (see Sec¬ 
tion 5.2.5). This is higher than what can be expected, since 
with reasonable assumptions that the explosion energy E e j < 
10£ bo i (Murase et al. 2014; Margutti et al. 2014) and a frac¬ 
tion £cr <0.1 of it going into CRs (Murase et al. 2011), the 
energy in neutrinos should be on the same order or less than 
£boi- Thus, also with a simple energetic argument, isotropic 
neutrino emission from PTF12csy causing the neutrino alert is 
implausible, especially on a time scale of seconds. However, 
a beamed emission from a jet with a small opening angle of 
< 30° would in principle be possible. 

4.2. X-ray Observations of PTF12csy 

The Swift satellite observed the supernova four times, on 
2012 April 20 (MJD 56037) and around 2012 November 15 
(MJD 56246) (see Table 2). We perform source detection us¬ 
ing the software developed for the 1SXPS catalog (Evans et al. 
2014) on each of the four observations, and on a summed im¬ 
age made by combining all the datasets. No counterpart to 
PTF12csy is detected. Upper limits are generated for each 
of these images, following Evans et al. (2014). A 28" radius 
circle centered on the optical position of PTF12csy is used 
to measure the detected X-ray counts c at this location, and 
the expected number of background counts Cb g , predicted by 
the background map created in the source detection process. 
We then use the Bayesian method of Kraft et al. (1991) to 
calculate the 3 cr upper limit cul on the X-ray count rate of 
PTF12csy, using the XRT exposure map to correct for any 
flux losses due to bad pixels on the XRT detector, and the 
finite size of the circular region. 

The upper limit count rate is converted to unabsorbed flux 
<I>ul using the HEASARC Tool WebPIMMS 4 , assuming a 
black body model with T = 0.6keV as in Miller et al. (2010b), 
a Galactic hydrogen column density of 1.31 x 10 21 cm~ 2 
(Willingale et al. 2013) 5 and a redshift of z = 0.0684. The re¬ 
sult is a 0.2 to lOkeV X-ray flux < 4.6 x 10 _14 ergcm _2 s“ 1 
for the most constraining upper limits, corresponding to a 0.2 
to lOkeV X-ray luminosity of Lx < 5.2 x 10 41 ergs 1 with 
a luminosity distance of about 308 Mpc. Using a power- 
law E 2 instead of a black body as an alternative X-ray 
emission model, the unabsorbed flux upper limit becomes < 

4 https://heasarc.gsf c.nasa.gov/Tools/w3pimms.html 

5 http://www.swift.ac.uk/analysis/nhtot/index.php 
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7.4 x 10 _14 ergcm _2 s“ 1 , and hence. Lx < 8.4 x 10 41 ergs" 1 . 

Comparing with other SNe Iln, e.g. SN 2008iy (Miller 
et al. 2010b) which had a measured X-ray luminosity of Lx = 
(2.4±0.8) x 10 41 ergs" 1 or SN 2010jl (Ofek et al. 2014b) 
with Lx ~ 1.5 x 10 41 ergs" 1 , we cannot exclude X-ray emis¬ 
sion from PTF12csy with our measured upper limit. How¬ 
ever, Svirski et al. (2012) suggest that Lx be about 10" 4 of 
the bolometric luminosity at the time of the shock breakout. 
With the estimated bolometric luminosity from Section 5.2.5 
around the time of the first Swift observations, this implies 
Lx ~ 6.4 x 10 38 ergs“ 1 , well below our X-ray limits. 

5. LOW-ENERGY FOLLOW-UP DATA 
5.1. Optical and UV Observations of PTF 12csy 

During the follow-up program of the neutrino alert, the first 
observations were done on 2012 April 03, 05, 07 and 09 (MJD 
56020 to 56026) by PTF with the Palomar Samuel Oschin 
48-inch telescope (P48) (Law et al. 2009), which is a wide- 
held Schmidt telescope. The images (see Fig. 2) revealed a so 
far undiscovered supernova, named PTFUcsy, at a magnitude 
of ~ 18.6 in the Mould //-band. More photometric observa¬ 
tions were carried out, with the P48 and the Palomar 60-inch 
(P60) telescopes (Law et al. 2009) at the Palomar Observa¬ 
tory in California, and the Faulkes Telescope North (FTN) 
(Brown et al. 2013) at Haleakala on Maui, Hawaii. Spec¬ 
troscopy was taken as well, with the Gemini North Multi- 
Object Spectrograph (GMOS) (Hook et al. 2004) on the 8 m 
Gemini North telescope (Mauna Kea, Hawaii) on 2012 April 
17 (MJD 56034) and with the Low-Resolution Imaging Spec¬ 
trometer (LRIS) (Oke et al. 1995) on the 10m Keck I tele¬ 
scope (Kamuela, Hawaii) on 2013 February 09 (MJD 56332), 
enabling the identification of the supernova as a Type Iln SN 
with narrow emission lines. The spectra are available from 
WISeREP 6 (Yaron & Gal-Yam 2012). 

P48 data were extracted using an aperture photometry 
pipeline and are calibrated with 21 close-by SDSS stars 
(Eisenstein et al. 2011; Gunn et al. 2006; Ahn et al. 2014; 
Alam et al. 2015). The faint host galaxy was subtracted and 
the upper limits are at the 5 cr level. P48 magnitudes are in 
the PTF natural AB magnitude system, which is similar, but 
not identical to the SDSS system. The difference is given by 
a color term, which is ignored in this work, except for the 
conversion of Mould R to SDSS r , explained in Section 5.2.1. 
The P60 photometry is tied to the same 21 SDSS calibration 
stars. Note that there might be host galaxy contamination in 
the late-epoch P60 photometry. P60 magnitudes are in the 
SDSS AB magnitude system. The FTN data were processed 
by an automatic pipeline, without host subtraction, and agree 
very well with the host-subtracted P60 data taken in the same 
night. 

The Pan-STARRS 1 (PS 1) telescope was not part of the real¬ 
time triggering and response system, but its wide-held cov¬ 
erage provides a useful archive to search retrospectively for 
detections. PS 1 is a 1.8 m telescope located at Haleakala on 
Maui in the Hawaiian islands, equipped with a 3.3° FOV and 
a 1.4 gigapixel camera (Kaiser 2004). In the course of its 3n 
steradian survey, the telescope observes each part of the sky 
typically 8 to 10 times per year (Magnier et al. 2013). PS1 
first detected PTF12csy on MJD 55847.582 and archived it 
as object PSO J104.6365+17.2622. The magnitudes in all 
PS1 images were obtained with PSF fitting within the Pan- 

6 http://wiserep.weizmann.ac.il 
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Figure 4. Background, gray, left axis: Gemini North spectrum from 2012 
April 17 (MJD 56034). Foreground, multiple colors, right axis: The filter 
response functions of the applied photometric filters, defined as filter trans¬ 
mission or effective area. Note that the absolute normalization is arbitrary 
and only the shape of the curves is relevant. 

STARRS Image Processing Pipeline (Magnier 2006). They 
are calibrated to typically seven local SDSS DR8 field stars. 
The magnitudes are in the natural PS 1 AB system as defined 
in Tonry et al. (2012), which is similar, but not exactly the 
same as SDSS AB magnitudes. Particularly the g-band can 
differ. 

The Swift UVOT data were analyzed using the publicly 
available Swift analysis tools (Nasa High Energy Astro¬ 
physics Science Archive Research Center (Heasarc) 2014) 7 , 
and a source was seen close to the detection threshold. 

ROTSE’s limiting magnitude of about 16 to 17 mag pre¬ 
vented a detection of the SN in ROTSE follow-up observa¬ 
tions. 

5.2. Photometry 
5.2.1. Photometric Corrections 

The photometry is corrected for Galactic extinction using 
R v =A V /E(B -V) = 3.1 and E(B-V) = 0.071 (Schlegel 
et al. 1998) 8 . The extinction coefficient is converted to the fil¬ 
ters’ effective wavelengths using the algorithm from Cardelli 
et al. (1989) 9 , i.e. 0.35mag for u, 0.30mag for B, 0.26mag 
for g, 0.19 mag for r, 0.11 mag for z. The extinction within 
the host galaxy could not be determined. 

In Fig. 4, the Gemini North spectrum is overlaid with the 
applied photometric filters. The strong Balmer lines con¬ 
tribute differently to the various filters. For the SED construc¬ 
tion (see Section 5.2.5), in order to approximate the black 
body continuum, the contribution of the strongest emission 
lines, Ha and H/L is removed from the photometry using 
the Gemini North spectrum and the filter curves. For Figs. 6 
and 7, the P48 Mould R magnitudes are converted to SDSS 
r by subtracting the Ha contribution (as above), applying the 
formulae in Ofek et al. (2012) valid for black body spectra, 
and then re-adding the Ha contribution to the r-band. After 
conversion, the P48 R magnitudes are consistent with the P60 
SDSS r magnitudes. 

The Swift UVOT data contain host contamination. Since 
no GALEX data from a pre- or post-SN epoch are available 

7 See http: //www. swift. ac.uk/analysis/uvot/ for instructions 

8 Obtained via the NASA/IPAC Infrared Science Archive http: //irsa. 
ipac.caltech.edu/applications/DUST/. 

9 Via http://dogwood.physics.mcmaster.ca/Acurve.html. 
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Figure 5. PTF12csy photometry in apparent magnitudes without applying 
corrections. The photometry is averaged over intervals of 10 days. The data 
originate from the following telescopes: uvw2, uvm2, uvw 1, u, b : UVOT; B: 
P60; g: P60, PS1, FTN; r : P60, PS1, FTN; R : P48; i : P60, PS1, FTN; z: P60, 
PS1; y: PS1. 

for the host galaxy 10 , no host subtraction can be done in the 
UV filters of UVOT. For the u, b and v filters, the host is sub¬ 
tracted by interpolating the host magnitudes from the SDSS 
DR12 data (Alam et al. 2015) 11 to the effective wavelengths 
of the UVOT filters. 

5.2.2. The Light Cun’es 

The earliest detection of PTF12csy was in the Pan- 
STARRS1 y-band on 2011 October 13 (MJD 55847.582), 169 
days prior to the neutrino alert in observer frame, correspond¬ 
ing to 158 days in host galaxy rest frame using z = 0.0684 
(see Section 5.3). The latest non-detection, again in Pan- 
STARRS1, was on 2011 March 21 (MJD 55641.3) in a 30 s 
Z-band frame, 206 days before the first detection (193 days in 
rest frame). Hence, the explosion time is not well constrained 
and can be anytime between MJD 55641.3 and MJD 55 847.6. 
Hereafter, we refer to the y-band detection at MJD 55 847.582 
as the first detection and use it as day 0 for the light curve. 

The uncorrected SN light curves with the data available 
through the IceCube optical follow-up program are displayed 
in Fig. 5, including photometry acquired with the Swift 
UVOT filters uvw2, mm2, uvwl, u and b\ the Johnson B filter 
on P60; the SDSS filters g, r, i with data from P60, PS1, and 
FTN; the SDSS z filter on P60; Mould R filter on P48; and 
Pan-STARRS y filter on PS 1. The entire uncorrected photom¬ 
etry in apparent magnitudes, as seen in Fig. 5, is also available 
in Table 3. The light curves are averaged within bins of 10 
days width, for each filter and telescope separately. Note that, 
in contrast to most of the other photometry, no host subtrac¬ 
tion is performed for the Swift UVOT magnitudes presented 
in Fig. 5 and Table 3. 

Figure 6 shows the light curve of selected filters in absolute 
magnitudes, after the photometric corrections. Light curves 
of other exceptional Type II SNe are overlaid for comparison: 
SN Iln 2006gy (Ofek et al. 2007; Smith et al. 2007; Kawabata 
et al. 2009; Miller et al. 2010a), one of the most luminous SNe 
ever recorded, and SN 2010jl (Stoll et al. 2011; Zhang et al. 
2012), a SN Iln that is spectroscopically similar to PTF12csy 

10 See http://galex.stsci.edu/GalexView/ 

11 http://skyserver.sdss.org/drl2/en/tools/chart/navi. 
aspx 



Figure 6. PTF12csy photometry (symbols) in absolute magnitudes, with cor¬ 
rection for Galactic extinction, and conversion of P48 Mould R magnitudes 
to SDSS r magnitudes (see Section 5.2.1). The data originate from the fol¬ 
lowing telescopes: g: P48, P60, PS1, FTN; r: P48, P6(XPS1, FTN; i: P60, 
PS1, FTN; z: P60, PS1; y: PS1. The photometry is averaged over intervals of 
10 days. Other absolute SN II light curves (lines) and a theoretical light curve 
from radioactive decay of nickel (black dashed line) are added for compari¬ 
son. The comparison light curves are partly not extinction corrected and have 
different reference dates (see text). 

(see Section 5.3) and shows signs of a collisionless shock in 
an optically thick circumstellar medium (CSM), hinting to¬ 
wards potential high-energy neutrino production (Ofek et al. 
2013). Note that the SN 201 Oj 1 light curve is not extinction 
corrected and the comparison light curves have different ref¬ 
erence dates: SN 201 Oj 1 is relative to maximum light, 2006gy 
is relative to the explosion time. A theoretical light curve from 
pure radioactive decay of 56 Ni —> 56 Co —> 56 Fe (black dashed 
line) is added to the figure as well, scaled to match the ob¬ 
served absolute magnitude of PTF12csy. 

The brightest observed absolute magnitudes after appli¬ 
cation of photometric corrections (see Section 5.2.1) and 
converting to absolute magnitudes with a distance modulus 
of fi = 37.443 (z = 0.0684) are M ? ss — 19.0mag, M r « 

— 19.0mag, Mi ss —19.6mag, M z « — 19.4mag, and M y ta 

— 19.0 mag, assuming standard cosmology with Hubble pa¬ 
rameter Hq = 70kms" 1 Mpc 1 , matter density Q. m =0.3, and 
dark energy density = 0.7. While these are lower lim¬ 
its to the peak magnitude due to the sparse sampling, these 
absolute magnitudes are relatively modest compared to the 
most luminous SNe Iln, e.g. SN 2006gy (Mr = —22 mag) 
(Kawabata et al. 2009) or SN 2008fz (My = —22.3 mag) 
(Drake et al. 2010). They are however comparable to the SNe 
Iln 2008iy (M r « -19.1 mag) (Miller et al. 2010b), 1988Z 
(Mr < —18.9mag) (Turatto et al. 1993) and SN 2010jl (Mr < 
—20.0mag) (Zhang et al. 2012). 

5.2.3. Decline Rates and Energy Source 

The light curves of PTF12csy indicate a plateau within 
~ 100days after first detection, and a slow fading afterwards. 
The corrected absolute magnitude light curves are fitted to ob¬ 
tain the linear decline rates in different photometric filters, 
during different epochs (see Fig. 7 and Table 4). For some 
epochs and filters, especially g and r, the decline rates are 
close to 0.98mag(100 days) -1 , the decline rate expected for 
radioactive 56 Co decay (Miller et al. 2010b), while in general 
decline rates are slower, indicating that at least part of the ra¬ 
diated energy is powered by interaction of the SN ejecta with 
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Table 3 

Photometric observations of PTF12csy 


MJD Date 

Rest frame days 

Mag 

Abs. mag 

Lim. mag 

Filter 

Tel. 

55 273.219 

-537.589 



21.11 

g 

P48 

55 294.162 

-517.987 



20.62 

R 

P48 

55431.514 

-389.429 



18.99 

R 

P48 

55477.402 

-346.478 



20.67 

R 

P48 

55 596.889 

-234.642 



18.90 

g 

P48 

55 641.304 

-193.071 



21.40 

z 

PS1 

55 847.588 

0.005 

18.52±0.08 

-18.92 


y 

PS1 

55 875.515 

26.145 

18.00±0.02 

-19.45 


i 

PS1 

55937.502 

84.163 

18.17 ± 0.04 

-19.27 


z 

PS1 

55948.816 

94.752 

18.82±0.05 

-18.62 


g 

PS1 

55948.841 

94.776 

18.66±0.02 

-18.79 


r 

PS1 

55957.475 

102.858 

19.00±0.02 

-18.44 


g 

PS1 

55967.269 

112.024 

18.23±0.01 

-19.22 


i 

PS1 

55997.366 

140.194 

18.80±0.08 

-18.64 


y 

PS1 

56022.981 

164.169 

18.61 ±0.08 

-18.84 

18.23 

R 

P48 

56026.246 

167.225 

19.01 ±0.08 

-18.43 


z 

PS1 

56034.841 

175.270 

19.47 ±0.04 

-17.98 


r 

P60 

56034.844 

175.273 

18.99±0.08 

-18.46 


z 

P60 

56035.177 

175.584 

18.64±0.04 

-18.80 


R 

P48 

56035.922 

176.281 

18.68 ±0.03 

-18.76 


i 

P60 

56035.925 

176.284 

20.15 ±0.06 

-17.30 


B 

P60 

56035.928 

176.287 

19.69±0.05 

-17.75 


g 

P60 

56036.581 

176.899 

18.66±0.04 

-18.78 

20.97 

i 

FTN 

56036.585 

176.902 

19.54±0.06 

-17.90 

21.12 

r 

ftn 

56036.590 

176.906 

19.73 ±0.07 

-17.72 

21.23 

g 

FTN 

56037.150 

177.431 



19.13 

V 

UVOT 

56037.150 

177.431 

20.60 ±0.35 

-16.84 

20.79 

u 

UVOT 

56037.150 

177.431 

22.21 ±0.28 

-15.23 

22.76 

uvm2 

UVOT 

56037.150 

177.431 

22.52±0.36 

-14.92 

22.71 

uvw2 

UVOT 

56037.150 

177.431 



21.84 

uvw 1 

UVOT 

56037.150 

177.431 

19.46±0.27 

-17.98 

19.95 

b 

UVOT 

56039.175 

179.326 

18.74±0.03 

-18.70 


i 

P60 

56039.177 

179.328 

19.55 ±0.04 

-17.89 


r 

P60 

56039.178 

179.329 

20.17 ±0.06 

-17.28 


B 

P60 

56039.181 

179.332 

19.73 ±0.04 

-17.71 


g 

P60 

56040.285 

180.365 

18.60±0.07 

-18.85 

20.98 

i 

FTN 

56043.178 

183.073 

18.70±0.05 

-18.74 


R 

P48 

56158.504 

291.015 

19.40 ±0.12 

-18.05 


i 

P60 

56161.496 

293.815 

20.58±0.12 

-16.86 


r 

P60 

56163.490 

295.682 

19.69 ±0.13 

-17.76 


i 

P60 

56165.486 

297.550 

21.01 ±0.15 

-16.44 


g 

P60 

56176.465 

307.826 

20.73 ±0.26 

-16.72 


g 

P60 

56177.453 

308.751 

20.07 ±0.26 

-17.37 


r 

P60 

56185.431 

316.218 

19.72 ±0.11 

-17.73 


i 

P60 

56185.432 

316.219 

20.77 ± 0.17 

-16.68 


r 

P60 

56185.436 

316.223 

20.92±0.15 

-16.53 


g 

P60 

56200.887 

330.685 

19.81 ±0.30 

-17.63 


i 

P60 

56202.385 

332.086 

21.09±0.25 

-16.35 


r 

P60 

56215.381 

344.250 

21.23±0.39 

-16.21 


r 

P60 

56215.429 

344.295 

20.11 ±0.10 

-17.33 


i 

P60 

56215.435 

344.301 

21.25±0.17 

-16.19 


g 

P60 

56219.344 

347.960 

21.16 ± 0.13 

-16.28 


g 

P60 

56219.837 

348.421 

20.97±0.13 

-16.47 


r 

P60 

56224.988 

353.242 

20.01 ±0.17 

-17.44 


i 

P60 

56229.808 

357.754 

20.20 ±0.45 

-17.24 


i 

P60 

56246.320 

373.208 

22.71 ±0.26 

-14.73 

23.34 

uvw2 

UVOT 

56246.320 

373.208 

22.55 ±0.35 

-14.89 

22.79 

uvw 1 

UVOT 

56246.320 

373.208 



23.05 

uvm2 

UVOT 


Note. — Rest frame days relative to first detection on MJD 55 847.582. 10-day binning. No correction for extinction. 


a dense CSM (Miller et al. 2010b). 

Additionally, radioactive decay of 56 Co at a still relatively 
high absolute magnitude of ~ —19 mag implies a preceding 
56 Ni decay with an extremely bright peak, which was not ob¬ 
served, although the data are quite sparse. Assuming that the 
luminosity is generated by radioactive decay alone and fol¬ 
lowing Kulkami (2005), we estimate that >1.7 M e of 56 Ni 
would be required to provide the bolometric luminosity of 
9.7 x 10 42 ergs _1 at lOOd in rest frame (see Section 5.2.5). 


The lower limit on the 56 Ni mass is set by assuming that the 
explosion and thus generation of 56 Ni was at the latest pos¬ 
sible time, directly before the first detection. Figure 6 shows 
the corresponding theoretical light curve resulting from the 
radioactive decay of nickel and cobalt (black dashed line). 
Adopting an earlier explosion time results in an even larger 
56 Ni mass. 

This is much more than the usual amount of 56 Ni of 
<O.5M 0 , often <0.1M S (see e.g. Pejcha & Thompson 
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Figure 7. Light curves of several filters with the fitted linear declines. See 
Table 4 for the numerical values of the found decline rates. 


Table 4 

Decline rates of the PTF12csy light curve 


Filter 

0to 150d a 

70 to 200 d a 

170 to 400 d a 

uvw2 



0.097 ±0.227 

g 


1.127 ±0.044 

0.893 ±0.051 

r 


0.974 ±0.053 

0.907 ±0.059 

i 

0.269 ±0.024 

0.656 ±0.089 

0.764 ±0.052 

z 


0.943 ±0.079 


y 

0.199 ±0.080 




a Units: mag(100d) 1 . Indicated periods in rest frame days 
relative to first detection on MJD 55847.582. 


(2015), also Margutti et al. (2014)). However, extremely su- 
perluminous SNe might have 56 Ni masses of that order of 
magnitude (Gal-Yam 2012). 

We note that in addition to the 56 Ni mass and luminosity 
arguments, the spectrum showing intermediate width Balmer 
lines and a continuum appears inconsistent with radioactive 
decay as well (see Section 5.3). 


steeper than —0.3 in order to reasonably describe the data. 
It lies in the range of —3 to —1.2, using the constraint on 
the explosion time (see Section 5.2.2) for the temporal zero 
point of the power-law. The best fit is at a = —3, with the 
explosion time at the lowest allowed value which is the date 
of the last non-detection. This suggests a very steep CSM 
density profile r~ 5 (see Ofek et al. 2014b, eq. 12), com¬ 
pared to the profile °e r~ 2 resulting from a wind with steady 
mass loss. However, the self-similar solutions of the hydro- 
dynamical equations (Chevalier 1982) that are used in Ofek 
et al. (2014b, eq. 12) are invalid if the CSM density profile is 
steeper than r“ 3 . But nevertheless, as discussed for the late¬ 
time light curve in Ofek et al. (2014b, sec. 5.2), probably the 
profile is steeper than r“ 3 . 

This leaves us with several possible explanations: 

1. Already between rest frame days 93 and 200, the SN 
was in the late, e.g. snow-plow, phase. This is consis¬ 
tent with SN 2010jl, where the late-time light curve also 
shows a power-law index a ss —3 (Ofek et al. 2014b, 
sec. 5.2). Assuming that the break in the light curve 
between power-law phase and late phase occurred just 
before the first r-band detection at day 93, and compar¬ 
ing with SN 201 Oj 1 (Ofek et al. 2014b), then this means 
that the SN was likely already a few hundred days old, 
and the r-band maximum was about 1 to 1.25 mag 
brighter than the observed one (cf. Ofek et al. 2014b, 
Fig. 1), consistent with SN 2010jl’s r-band maximum. 
It follows that the power-law phase ended < 286 rest 
frame days after explosion, from which we can derive a 
swept CSM mass of < 12M 0 , using Ofek et al. (2014b, 
eq. 22) and adopting the standard values given for SN 
2010jl. 

2. The SN is powered by ejecta-CSM interaction, but its 
light curve is declining steeper than a f~° 3 power-law. 
This is possible, e.g. if spherical symmetry, assumed in 
Ofek et al. (2014b), is broken, if the optical depth is 
lower than in SN 201 Oj 1 or if the CSM density profile 
falls steeper than r~ 2 (s.a.). 


5.2.4. Fitting to an Interaction Model 

Here we assume that the light curve is powered by con¬ 
version of the ejecta’s kinetic energy to luminosity through 
interaction of the ejecta with the CSM. Following Ofek et al. 
(2014b, see also: e.g. Chugai & Danziger (1994); Svirski 
et al. (2012); Moriya et al. (2013)), we model the light curve 
as a power-law of the form L(t) = Ltf a . 

After shock breakout, there is a phase of power-law decline 
of the luminosity, with an index of typically a « —0.3. This 
lasts until the shock runs over a CSM mass equivalent to the 
ejecta mass and the shock enters a new phase of either conser¬ 
vation of energy if the density is low enough and the gas can¬ 
not cool quickly (the Sedov-Taylor phase), or conservation of 
momentum if the gas radiates its energy via fast cooling (the 
snow-plow phase). During the late stage, the light curve will 
be declining more steeply, in both cases (Ofek et al. 2014b). 

Since we assume PTF12csy to be powered by interaction, 
we try to fit the interaction model from Ofek et al. (2014b) 
to the light curve data with the least-squares method. We per¬ 
form the fit within the range of 93 to 200 rest frame days, start¬ 
ing at the first r-band detection, and use the r-band light curve 
scaled with the bolometric luminosity from Section 5.2.5. It 
is found that the power-law index a needs to be significantly 


3. The SN is not powered by interaction, but by radioac¬ 
tive decay, leading to an exponential light curve decline. 
However, this appears unlikely, as noted above in Sec¬ 
tion 5.2.3. 

5.2.5. Spectral Energy Distribution (SED) 

Since the spectra are only roughly calibrated, the spectral 
energy distribution (SED) is approximated from photometric 
data. For the highest spectral range and number of observa¬ 
tions, a window of 10 observer frame days around day 189, 
from day 184 to 194, is used to select data (day 172.2 to day 
181.6 in rest frame). Photometric corrections are applied, e.g. 
Ha and H/j removed (see Section 5.2.1). The data are plot¬ 
ted in Fig. 8 as function of the filters’ effective wavelengths. 
A black body spectrum is assumed to describe the SED and 
is fitted to the data. For each filter, a model data point cor¬ 
responding to the black body spectrum is calculated via in¬ 
tegration of the black body spectrum and the filter function 
following the SDSS definition of AB magnitude in Fukugita 
et al. (1996, eq. 7). A "/} fit minimizes the difference between 
the model data and the measured data. 

The fit results in a reduced X 1 / n Aof = 7.9/5 = 1.6 and de¬ 
livers estimates for both the rest frame temperature T and the 
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Observed Wavelength [A] 


Figure 8. SED of PTF12csy using photometry from 10 days around 
day 189 (observer frame) after the first detection. The fitted rest frame 
temperature is T = (7160±270)K and the fitted bolometric luminosity 
(5.53± 1.18) x 10 42 ergs' 1 . 

absolute bolometric luminosity Lbol of the photosphere emit¬ 
ting the black body radiation: T = (7160 ± 270) K and L| 10 i = 
(5.53 ± 1.18) x 10 42 ergs ~ 1 , where the errors correspond to 
1 a. Applying the Stefan-Boltzmann law, we can calculate 
the radius of the black body photosphere from the bolometric 
luminosity. We estimate it to be R v hot = (1-7 ±0.1) x 10 15 cm. 

Finally, to obtain an estimate on the total radiated energy, 
the lines’ contributions to the luminosity have to be added to 
the continuum luminosity. Using the Gemini North spectrum, 
the contribution of the Ha and H/i line to the total luminosity 
is computed and added to the continuum luminosity from the 
black body fit. This results in an estimated total radiated lumi¬ 
nosity of (6.4± 1.2) x 10 42 ergs _I at day 189 in the observer 
frame, i.e. day 177 in the rest frame. 

We use the fitted shape of the /-band light curve (cf. Sec¬ 
tion 5.2, Fig. 7, Table 4) to extrapolate this value and find 
a total radiated luminosity of ~ 9.7 x 10 42 erg s 1 at lOOd 
(rest frame), as used in Section 5.2.3, and a total energy 
of T’hoi = 2.1 x 10 50 erg radiated within 400 rest frame days 
after first detection, comparable to SN 2008iy which had 
~ 2 x 10 50 erg) (Miller et al. 2010b) and SN 2010jl with 
4.3 x 10 50 erg (Zhang et al. 2012). This is a lower limit on the 
total radiated energy, since we lack photometric data between 
explosion and first detection and do not extrapolate before the 
first detection. Additionally, as discussed below, we are ne¬ 
glecting a possible contribution of X-ray and y-ray emission 
to the total radiated energy, which is not considered by the 
black body spectrum based on the UV and optical data. 

We recommend to treat these results with caution, since 
Ofek et al. (2014b) pointed out that at late times the fraction of 
energy released from SNe Iln in X-rays can increase, causing 
the optical spectrum to deviate from a black body as fewer 
photons are available in the optical. This can lead to an ef¬ 
fective decrease of the estimated photospheric radius. In this 
context, our estimates of /? p hot- Aral. and Aral must be treated 
as lower limits. Unfortunately, the X-ray flux from PTF12csy 
was not detected (see Section 4.2). 

5.3. Spectroscopy 

Two spectra were acquired (Table 5, Fig. 9). They are 
dominated by narrow emission lines, characteristic for Type 
Iln SNe, with a very weak blue continuum, which indicates 


Table 5 

Log of spectral observations 


MJD 

Tlisc 

Tdet 

Av(kms ’) 

Instrument 

56034 

ti 

175 

80 

Gemini North GMOS 

56332 

290 

454 

100 

Keck I LRIS 


Note. — 7d isc : Rest frame days after PTF discovery. 7d e u Rest frame days 
after first detection by PS1. Av: spectral resolution at the Ha line at 7014 A 
in observer frame. 


the old age of the SN. No continuum is visible in the late 
spectrum. The SN emission lines are primarily hydrogen, 
the Baimer series is visible from Ha up to He. The oxygen 
lines OI AA7772,7774,7775,8447, Oil A3727 with FWHM 
« 500kms-\ and Om AA4364,4960,5008 with FWHM 
« 350km s are very narrow and were most likely produced 
by circumstellar gas released by the progenitor prior to ex¬ 
plosion and then photoionized by UV radiation (Filippenko 
1997). 

Figure 10 shows a close-up on the Ha line from both spec¬ 
tra, plotted vs. Doppler velocity relative to the rest frame 
line center. In the early spectrum, the Ha line peaks at the 
line center and the line is composed of a narrow, intermedi¬ 
ate and broad component with FWHM of ~ 400 km s *, ~ 
200010ns" 1 , and ~ 5000kms 1 respectively, found by fitting 
a superposition of three Gaussian functions to the Ha profile. 
This is similar to other SNe Iln, e.g. SN 1988Z and 2008iy 
(Turatto et al. 1993; Filippenko 1997; Miller et al. 2010b). 
The Ha profile with broad and intermediate-width component 
can be explained as a result of the interaction of the SN ejecta 
with a two-component wind (Chugai & Danziger 1994). In 
this model, the broad line component is emitted from shocked 
SN ejecta expanding in a relatively rarefied wind, while the 
intermediate component arises from a shocked dense part of 
the wind, which can either consist of dense clumps or be a 
dense equatorial wind (Chugai & Danziger 1994). This is an 
indication of ejecta-CSM interaction. 

The intermediate and broad component of the early spec¬ 
trum’s Ha line are blueshifted which may indicate formation 
of dust, as explored by Smith et al. (2012) for SN 2010jl. A 
more recent study by Gall et al. (2014), the most comprehen¬ 
sive work on the SN 2010jl emission line blueshifts to date, 
finds very strong evidence for a wavelength dependence of 
the blueshift. Therefore, the authors conclude that the origin 
of the blueshifts is most likely the rapid formation of large 
dust grains, confirming Smith et al. (2012) and having impli¬ 
cations on the origin of dust in galaxies. 

Alternatively, Fransson et al. (2014) explain the line 
blueshift in SN 201 Oj 1 with a bulk velocity of the emitting 
gas towards the observer. This is more consistent with obser¬ 
vations if the spectral lines are symmetric about a center and if 
there is no wavelength dependence of the blueshift. The bulk 
velocity is believed to be the result of radiative acceleration 
of the gas by flux from the SN. Presumably, there are also 
other possible explanations for the blueshift, e.g. the geome¬ 
try or density structure of the CSM. In case of SN PTF12csy, 
spectral line blueshift is only clearly visible in the Ha line, 
prohibiting the interpretation of the blueshift in favor of any 
scenario. 

The late spectrum’s Ha line has a much more complicated 
structure than the early spectrum’s. It does not peak at zero 
velocity anymore, but the peak is blueshifted and there are 
many sub-peaks. Again, the blueshift might be connected to 
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Figure 9. Spectra taken with Gemini North on 2012 April 17 (top) and Keck I on 2013 February 09, showing narrow (Type Iln) emission lines. The Ha line at 
~ 700()A (observer frame) is the strongest emission line and has a complicated structure. See Fig. 10 for a close-up of the Ha line. 



Figure 10. Comparison of the Ha line in both spectra. The x-axis shows the 
Doppler velocity relative to the line center at 6564.61 A, assuming a redshift 
of 0.0684. 

dust formation or radiative gas acceleration, but other reasons 
are conceivable as well. At least part of the late Ha line’s 
complex appearance might be due to superposition of other 
spectral lines, e.g. N II. Apart from that, it is an indication of 
an inhomogeneous, maybe clumpy, CSM structure, and per¬ 
haps asymmetric SN explosion. 

The spectra are compared to template spectra from the 
Padova-Asiago Supernova Archive (ASA) (Harutyunyan et al. 


2008) using the online tool GELATO 12 . The algorithm (Haru¬ 
tyunyan et al. 2008) divides a spectrum into 11 relevant bins 
and averages within the bins to classify and compare with 
the archived spectra. The PTF12csy spectra are de-reddened 
with E(B — V) =0.1 and compared to other Type II SNe. 
GELATO returns the best 30 matching spectra together with 
their phases, ordered by quality of fit. For both the Gemini 
North spectrum taken at 175 d and the Keck I spectrum from 
454 d, the majority of best matching template spectra come 
from SN 2010jl. The mean phase of the matching spectra is 
significantly higher for the Keck I spectrum: (378 ± 102) d, 
versus (154 ± 21) d for the Gemini N spectrum. However, the 
reference dates for the spectra are mostly the discovery date, 
only rarely the date of maximum light or explosion date. 

5.4. Host Galaxy 

The galaxy hosting PTF12csy is a faint dwarf galaxy des¬ 
ignated SDSS J065832.82+171541.6 13 , barely visible in the 
SDSS DR12 images. Since it has no cataloged redshift, we 
measure its redshift via Oil A3727 and Om AA4960,5008 
lines in the Keck I spectrum. The redshift is 0.0684 ± 0.0001, 
corresponding to a luminosity distance of ~ 300Mpc as¬ 
suming standard cosmology with Hubble parameter Hq = 
70kms -1 Mpc -1 , matter density Q. m = 0.3, and dark energy 
density Q A = 0.7. 

12 https://gelato.tng.iac.es 

13 http://skyserver.sdss.org/dr12/en/tools/explore/ 
summary.aspx?id=0xll2dlf06c01f0a28 
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Adopting the luminosity distance from above, the host 
galaxy has absolute magnitudes of M g —16.2mag, M r ta 
— 16.6 mag and M, « —16.7 mag 13 (corrected for Galactic ex¬ 
tinction). Using the luminosity-metallicity relation from Lee 
et al. (2006, eq. 1), we find a metallicity of 12 + logO/H « 8, 
indicating that the host galaxy is quite metal-poor. Over- 
luminous SNe Iln, such as PTF12csy, have been preferen¬ 
tially found to occur in subluminous, low-metallicity galax¬ 
ies (Miller et al. 2010b; Stoll et al. 201 1), such as the host of 
PTF12csy. This is a trend, which is also observed for long 
GRBs (Stoll et al. 2011). Statistics are still low and Miller 
et al. (2010b) cautioned that there could be some selection 
bias due to intrinsically bright SNe in faint host galaxies be¬ 
ing more easily discovered during surveys doing aperture pho¬ 
tometry. However, new surveys performing image subtrac¬ 
tion and observing large untargeted fields, e.g. PTF and Pan- 
STARRS, provide increasing evidence for this trend, as most 
of the discovered bright objects would have been luminous 
enough to be detected in bright galaxies and in searches that 
are targeted to bright galaxies (Stoll et al. 2011). PTF12csy, 
probably discovered by coincidence in an unbiased way, con¬ 
firms this emerging trend as well. 

The SDSS DR12 1 has the host galaxy’s cataloged center 
position about 3" away from the found SN position (see also 
Fig. 2). With an apparent diameter of about 5 ", corresponding 
to ~ 7kpc, this is quite far from the center of the galaxy, i.e. 
about 4kpc off-center. 

6. SUMMARY AND CONCLUSION 

The highest significance alert from the IceCube Optical 
Follow-Up program led to a coincidental discovery of the in¬ 
teresting and unusual Type Iln SN PTF12csy, which was al¬ 
ready at least 169 days old. The combined a posteriori sig¬ 
nificance of the neutrino doublet alert and the coincident de¬ 
tection of any core-collapse SN within the error radius of the 
neutrino events (0.54°) and within the luminosity distance of 
the SN (300Mpc) is 2.2 a, for the time interval of the IceCube 
data acquisition season 2011/12. 

PTF12csy is rare and unusual: With peak absolute magni¬ 
tudes of M r < —19, perhaps about -20, it belongs to the most 
luminous SNe. The SN is most likely powered by interac¬ 
tion of the ejecta with a dense circumstellar medium (CSM). 
The spectrum indicates a complicated structure of the CSM. 
Its host galaxy is a faint and metal-poor dwarf galaxy, con¬ 
firming an observed trend for luminous SNe Iln. PTF12csy 
is similar in photometry and spectroscopy to other rare lumi¬ 
nous SNe Iln, e.g. SNe 2008iy and 2010jl. The total radiated 
energy is 2 x 10 50 erg within the first 400 rest frame days after 
detection. 

Given the ejecta-CSM interaction, high-energy (HE) cos¬ 
mic ray production and neutrino emission may be expected 
on a time scale of 1 to lOmonths, according to Murase et al. 
(2011) and Katz et al. (2012). However, the SN is too far away 
for IceCube to detect this emission. A complementary neu¬ 
trino analysis performed offline, using one year of IceCube 
data which cover most of the optical SN fluence, does not re¬ 
veal a signal-like accumulation of neutrino events from the 
SN’s position, leading to a very high upper limit of more than 
1000 times the tested model fluence, owing to the large dis¬ 
tance. 

Due to the long delay of several months between explo¬ 
sion date and neutrinos, the doublet of neutrinos within less 
than two seconds cannot be explained by the formation of 


a jet shortly after core-collapse according to the choked jet 
model (Ando & Beacom 2005). Nor can it be explained by 
the expected HE neutrino production from ejecta-CSM in¬ 
teraction of SNe Iln (s.a.). We therefore conclude that the 
only reasonable explanation is that the SN detection was co¬ 
incidental and the neutrino doublet was produced by uncorre¬ 
lated background events of atmospheric neutrinos and/or mis- 
reconstructed atmospheric muons 14 . 

However, if there was a delaying mechanism at work, such 
as the spin-down of a supramassive neutron star that delays 
its collapse to a black hole, as in the blitzar model (Falcke & 
Rezzolla 2014), then the neutrino doublet may have originated 
from a jet forming in the course of this delayed collapse few 
hundred days after SN explosion. 

The coincidental detection of a Type Iln SN following 
an IceCube neutrino alert demonstrates the capability of the 
follow-up system to reveal transient HE neutrino sources. An 
advantage of the follow-up paradigm is the prompt availabil¬ 
ity of multi-messenger information for the identification of 
the source, as well as the mere statistical significance of a 
coincidence between a neutrino burst and an electromagnetic 
transient detection. In this case, the significance is very low 
due to the delay of several months between explosion and neu¬ 
trinos. However, this coincidence motivates the continuation 
of the follow-up program as well as further stacked neutrino 
analyses of Type Iln supernovae. 
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